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摘 要 : 降雨 波动 是 荒漠 植被 生长 发 育 最 主要 制约 因子 , 直 


SAA 


ES Ml 3) AE BALM 结构 和 功能 的 变化 。 梭 梭 作为 干旱 


区 植被 的 优势 种 之 一 ,在 防风 固沙 涵养 水 源 方面 


力 对 降雨 事件 的 响应 过 程 。 本 文 以 黑河 中 游荡 漠 - RY 


扮演 重 


角色 ,本 文 的 目的 是 理解 不 同 林 龄 人 工 梭 梭 林 植被 生产 


川 过 渡 带 不 同 林 龄 人 工 梭 梭 为 研究 对 象 ,以 2017 一 2018 年 5 
1 日 至 9 月 31 日 时 序 MODIS-NDVI( 近 似 于 年 生产 力 ) 日 数据 和 对 应 日 降雨 数据 为 基础 ,采用 阅 值 模型 和 统计 学 


方法 分 析 了 其 NDVI 对 不 同 降雨 等 级 的 滞后 响应 过 程 。 结 果 表 明 :Q@ 降雨 小 于 2 mm 时 ,不 同 林 龄 梭 权 NDVI 增长 
率 均 在 5% 以 上 ,大 于 40 a 生 的 反应 最 迅速 ,增长 率 为 10% ~56% ;2 ~5 mm 时 ,10 ~20 a 生 的 变 幅 最 大 ,增长 率 为 
11% ~83% ;5 ~10 mm 时 ,20 ~40 a 的 最 为 敏感 ,增长 幅度 最 高 达到 170% ;大 于 20 mm 时 ,大 于 40 a 的 NDVI 增长 


率 大 于 其 他 3 个 林 龄 ,增幅 最 高 达 


76% . D 降雨 后 ,不 同 林 龄 人 工 梭 梭 NDVI 响应 时 间 不 同 ,0 ~ 10 a 的 为 7.4 


( £2.8)d,10 ~20 a 的 为 8.3( +3.1)d,20 ~40 a 的 为 8.1( +2.7)d, KF 40 a 的 为 8.2( +3.2)d。@ 降雨 脉动 是 
MERER - 绿洲 过 滤 带 植被 生产 力 迅速 变化 的 关键 因素 。 不 同 林 龄 梭 梭 对 不 同 级 别 降雨 响应 模式 不 同 ,2 ~ 10 mm 
降雨 条 件 下 ,人 工 梭 梭 林 NDVI 增 加 幅度 最 大 ,高 于 其 他 2 个 降雨 量 量 级 NDVI 响应 变化 幅度 。 依 据 上 述 分 析 可 以 


推断 ,在 蕊 并 -绿洲 过 渡 
的 作用 是 最 重要 的 。 
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在 全 球 气候 变化 背景 下 ,中 国 西部 干旱 区 主要 
以 气温 和 降水 量 增加 为 主要 响应 方式 " ,其 中 依赖 
冰雪 融 水 和 降水 维持 的 干旱 区 绿洲 和 荒漠 生态 系统 
是 响应 变化 最 敏感 的 生态 系统 " 。 特 别 是 对 于 荡 
漠 -绿洲 交错 区 的 旱 生 型 植被 ,由 于 水 分 来 源 主要 
依赖 降水 ,所 以 该 地 区 的 植被 水 分 生理 代谢 ,生长 发 
育 和 系统 生产 力 等 不 同 尺度 生态 过 程 都 受 生 长 期 降 
水 大 小 和 频率 的 影响 。 梭 梭 作为 干旱 区 最 主要 
的 国 沙 灌木 被 广泛 应 用 于 荒漠 - 绿洲 过 渡 带 的 人 工 
植被 防护 体系 建设 ,因此 ,研究 不 同 林 龄 梭 梭 对 降水 
响应 模式 是 对 干旱 区 人 工 梭 梭 林 水 分 适应 机 制 的 深 
入 了 解 ,也 是 干旱 区 生态 水 文学 研究 重要 方向 之 一 。 

有 关 干 旱 区 植被 与 水 分 关系 , Huang 等 中 在 西 
北 温带 幕 漠 生态 系统 通过 生长 季 典 型 灌木 林地 添加 
氨 实 验 , 评 估 了 降雨 对 该 生态 系统 净 交 换 量 (net e- 
影响 ,结果 发 现 
NEE 对 降雨 的 响应 程度 主要 依赖 于 降雨 时 间 和 植 
物 生长 节律 的 耦合 关系 。Liu 等 中 在 荒漠 生态 系统 


cosystem carbon exchange, NEE ) 的 
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带 以 10 mm 以 下 降雨 为 主导 的 条 件 下 ,10 ~20 a 生 及 20 ~40 a 生 人 工 梭 梭 林 在 防护 体系 


的 研究 表明 ,生长 季 前 期 足够 的 降水 导致 更 多 一 年 
生 浅 根植 物 的 萌发 ,生长 季 降 雨量 与 频率 都 显著 影 
响 生 态 系统 碳 储量 的 年 际 变化 。 从 与 梭 梭 相关 研究 
来 看 , 梭 梭 幼苗 在 不 同 梯度 供水 处 理 下 抗旱 性 存在 
明显 差异 ,供水 越 多 越 有 利于 梭 梭 幼苗 存活 与 生长 ， 
随 供水 量 增加 梭 梭 根 长 . 干 重 . 地 上 生物 量 均 迅 速 增 
KET, Liu 等 中 在 准噶尔 例 地 通过 同位 素 方法 测 
定 了 梭 梭 水 分 来 源 和 水 分 利用 策略 ,发 现 地 下 水 、 地 
表 融 雪 形 成 的 浅 层 土壤 水 及 生长 季 降 雨 是 梭 梭 水 分 
的 直接 来 源 , 在 单 次 降雨 后 4 d 对 植物 枝条 水 分 样 
本 进行 测定 显示 ,木质 部 的 S” O 值 急 剧 增加 ,表明 
梭 梭 对 降雨 利用 具有 滞后 性 ,降雨 是 梭 梭 的 水 分 来 
源 之 一 。 从 梭 梭 林地 水 分 利用 供给 关系 转换 来 看 ， 
在 一 些 沙 漠 和 绿洲 过 渡 地 区 ,由 于 地 下 水 过 度 开采 
导致 地 下 水 位 显著 下 降 , 以 梭 梭 为 主 的 固沙 植物 由 
主要 利用 地 下 水 生长 逐渐 转变 为 雨 养 植物 ”0 。 
这 些 研究 表明 , 梭 梭 水 分 利用 策略 在 不 同 生 育 期 .不 
同时 间 随 水 分 供给 关系 变化 不 断 进行 权衡 调整 ,其 
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中 对 有 效 降雨 的 依赖 是 其 主要 调整 策略 。 此 外 ,一 
些 研究 发 现在 空气 湿度 相对 丰沛 的 6 一 8 月 ,凝结 水 
对 植物 水 分 生理 代谢 有 一 定 影响 ,但 对 生产 力 的 影 
响 尚 不 明确 。 最 后 ,从 大 尺度 气候 变化 趋势 来 
看 ,在 过 去 50 a 中 亚 干 旱 区 降雨 呈 增 加 趋势 |。 
因此 ,可 以 推断 在 地 下 水 位 持续 下 降 、 凝 结 水 供给 意 
义 不 大 的 前 提 下 ,降水 量变 化 是 影响 不 同 林 龄 人 工 
(或 自然 ) 梭 梭 林 生长 发 育 的 关键 ,其 对 降水 的 响应 
模式 研究 是 人 工 权 梭 林 更 新 管理 与 荒漠 -绿洲 过 渡 
带 人 工 植被 稳定 性 维持 的 重要 保障 。 

据 此 ,本 文 以 2017 一 2018 年 5 月 1 日 至 9 月 31 
日 时 序 MODIS-NDVI 日 数据 为 基础 ,以 不 同 林 龄 人 
工 梭 梭 林 (0 ~10 a,10 ~20 a,20 ~40 a, 大 于 40 a) 
为 研究 对 象 ,分 析 不 同 林 龄 人 工 权 权 对 降雨 脉动 事 
件 的 响应 模式 ,以 期 为 荒漠 - 绿洲 过 渡 带 人 工 防 护 
体系 的 栽植 与 维护 提供 科学 依据 。 


1 材料 与 方法 


1.1 研究 区 概况 

研究 区 位 于 中 国 科 学 院 临 泽 内 陆 河 流域 综合 研 
究 站 东部 5 km 张掖 临 泽 高 台 绿 洲 与 巴 丹 吉林 沙漠 
HIRE -绿洲 过 渡 带 。 该 区 属 温带 大 陆 性 气候 ,年 
平均 气温 7.6 °C ,年 平均 降水 量 117 mm, 其 中 70% 
左右 集中 在 5 一 9 月 。 该 区 自然 植被 组 成 主要 有 油 
$ (Artemisia ordosica) 、 红 人 砂 ( Reaumuria songarica ) 、 
珍珠 柴 ( Salsola passerine ) | Mi JA GE ( Eragrostis pilo- 
sa) Z VK ( Bassia dasyphylla ) 和 泡 泡 刺 ( Nitraria 
sphaerocarpa ) 等 ,人工 固沙 植被 主要 是 梭 梭 林 
(Haloxylon ammodendron) 。 经 查阅 临 泽 站 日 志和 当 
年 参与 造林 人 员 现 场 确 认 , 本 研究 梭 梭 林 林 龄 划分 
4 组 ,分 别 为 小 于 10 a、10 ~20 a 20 ~40 a 和 大 于 40 
a( 图 1, 表 1)。 
1.2 数据 源 及 处 理 
1.2.1 NDVI 数据 ”本 研究 基于 NASA 数据 中 心 
(https://ladsweb. modaps. eosdis. nasa. gov/search/ 
order/1 ,轨道 号 为 h25v05 ) 从 可 见 光 到 短波 红外 前 7 
个 波段 的 250 m 地 表 反 射 率 数 据 , 选 取 2017 .2018 
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图 1 研究 区 位 置 示意 图 
Fig.1 Location of the study area 
表 1 不 同 林 龄 梭 权 形 态 特征 
Tab.1 Morphological characteristics of Haloxylon 
ammodendron in different stand ages 
林 龄 /a ”密度 高 度 /m 基 径 /cm 冠 幅 /m? 
2 2mx2m 0.59+0.04d 0.97+0.08d 0.08 +0.05e 
5 lmx2m 2.02+0.10c 5.10+0.37c 3.02 +0.35d 
9 lmx2m 3.19+0.15b 8.45+0.58b 6.78+1.1lc 


13 lmx2m 3.41 +0.17b 9.85+0.69ab 10.48 +2.41be 
16 1 mx2m 3.75+0.15b 9.93+0.79ab 11.31 +1.09abc 
31 1 mx2m 3.93+0.17ab 10.17+0.58ab 12.02 +0.96ab 
41 1 mx2m 4.78+0.24a 12.99+1.2la 20.20 +2.93a 


注 :数据 来 自 文献 [13] (引用 张 珂 的 文献 作为 梭 梭 林 植 被 背景 ) 。 


对 Red 和 NIR 波段 进行 线性 变换 (NIR - Red )/ 
(NIR + Red) ,获取 2017 .2018 年 5 月 1 日 至 9 月 30 
日 的 NDVI 日 数据 集 。 
1.2.2 气象 数据 日 降雨 数据 来 自 于 中 国 科学 院 
临 泽 内 陆 河 流域 综合 研究 站 ,2017 一 2018 年 5 一 9 
月 降雨 共 29 次 ,参照 李 芳 等 ”的 分 类 标准 ,将 降 十 
量 大 小 分 为 5 个 等 级 :0 ~2 mm、 2 ~5 mm,5 ~ 10 
mm,10 ~20 mm, >20 mm。 独 立 降雨 事件 选择 :经 
历 降雨 过 程 后 , NDVI 值 平均 增加 10% 以 上 则 视 为 
有 效 降 雨 。 
1.3 分析 方 法 

阅 值 延迟 模型 。 模 型 存在 的 假设 前 提 是 对 


年 5 月 1 日 至 9 月 30 日 的 MODIS 日 地 表 反 射 率 产 
品 (Terra: MODO9GQ ) ,数据 处 理 先 用 MODIS 投影 
转换 软件 MRT( MODIS Reprojection Tools) 进行 重 
采样 和 投影 变换 , 重 采样 分 辨 率 为 250 m ,最终 输 出 
为 WGS -1984 坐标 系 下 的 CeoTif 格式 文件 。 通 过 


降雨 产生 响应 的 量 级 有 上 限 ,该 模型 可 用 来 描述 植 
物 生长 对 降雨 等 的 延迟 响应 ,模型 表达 式 为 : 


y=k Yy, +6, (1) 


式 中 :y 是 在 上 时 的 NDVI 庆 表示 NDVI 变化 速率 ; 


= 
a 
= 
< 
=a 
` 
($) 
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y,,At-1 时 的 NDVI;6, 表示 NDVI 在 降雨 后 的 


增 量 ; 


6, = Min | Yn x (1 =k) ,6 x (1 =) | © 


式 中 :yi 是 啊 应 变化 中 的 最 大 值 ;6 为 在 1 时 的 最 
大 增 量 。 
One (R - R’) R’ <R < RY 
R! — R} t t-T 
ô, = 0 R,_.<R’ (3) 
Ô max > R,_, >R” 


P Sn ERK RTE N E; RU 和 尼 分 别 为 
降雨 量 的 上 限 值 和 下 限 值 ;7 是 某 次 降雨 事件 后 
NDVI 发 生 响 应 的 时 间 间 隔 ;R,_, 是 有 效 降雨 事件 。 
在 阔 值 模型 中 ,各 项 参数 由 降雨 事件 导致 的 ANPP 
大 于 10% 的 降雨 事件 求 算 。 

利用 SPSS22 进行 单 因 素 方差 分 析 (ANOVA ) ， 
比较 不 同等 级 降雨 条 件 对 不 同 林 龄 人 工 梭 梭 林 植 被 
生产 力 影 响 的 显著 性 。 


2 结果 与 分 析 


2.1 降雨 格局 

从 研究 期 降雨 量 分 布 格局 来 看 ,大 部 分 降雨 以 
<5 mm 为 主 , 约 占 整个 生长 生 (5 一 9 月 ) 内 降雨 事 
件 的 51.7% (图 2) ,平均 单 次 降雨 量 为 1 mm;5 ~ 10 
mm 降雨 约 占 20.7% ,平均 单 次 降雨 量 为 7.8 mm; 
KF 10 mm 的 降雨 发 生 频率 最 低 , 约 占 13. 8% , 平 
均 单 次 降雨 量 为 28.5 mm。 其 中 2018 年 8 月 31 日 
至 9 月 1 日 单 次 降雨 量 达 到 34.6 mm, 

2017 年 和 2018 年 5 一 9 月 累积 降雨 量 分 别 为 
105.6mm 和 86.5mm( 图 3) ,2017 年 生长 季 5 个 等 
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图 2 降雨 等 级 与 降雨 频率 的 关系 


Fig.2 Relationship between rainfall level and rainfall frequency 


级 降雨 所 占 比 例 分 别 为 27% . 18% 27% 、9% 和 
18% ;2018 年 的 单 次 降雨 量 大 多 较 小 ,小 于 2 mm 降 
雨 等 级 所 占 比例 为 66.7% ,但 降雨 量 只 占 生长 季 降 
雨量 的 14. 2% 。 生 长 季 总 降雨 量 中 ,大 部 分 来 自 
>10 mm 降雨 。 所 有 这 些 特征 都 符合 干旱 区 降雨 在 
时 间 和 等 级 上 的 分 布 特征 。 
2.2 AMBRE NDVI 对 降雨 的 响应 

从 阔 值 延迟 模型 分 析 结 果 来 看 ( 表 2) ,不 同 林 
龄 人 工 梭 梭 林 NDVI 变化 速率 (k 值 ) 分 别 为 0. 384 
( 林 龄 小 于 10 a) .0.836( 林 龄 10 ~20 a) 1.706( 林 
内 20 ~40 a) 和 0.892( 林 龄 大 于 40 a) , NDVI 的 变 
化 速率 和 最 大 增长 量 (5。.) 分 别 在 20 ~40 a 和 10 ~ 
20 a 达到 最 大 ,二 者 变异 系数 分 别 为 42. 5% 和 
39.7% ;而 不 同 林 龄 NDVI 的 最 大 响应 量 (y,,, ) 和 降 
雨 响应 时 间 (7) 相差 不 大 ,变异 系数 分 别 为 10. 6% 
和 5.1% 。 说 明 不 同 林 龄 梯度 对 模型 4 个 参数 中 的 
k 和 6 影响 较 大 ,而 yy 和 7 受 影响 较 小 。 

从 NDVI 对 降雨 事件 的 响应 速率 来 看 ,降雨 事 
件 后 不 同 降雨 等 级 和 不 同 林 龄 人 工 梭 梭 林 对 应 的 
NDVI 增 长 率 也 明显 不 同 。 当 降雨 量 小 于 2 mm 时 
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1 9 17 25 33 41 49 57 65 73 81 89 97 105 113 121 129 137 145 153 
天 数 /d ( 1 为 5 月 1 日 ) 
图 3 生长 季 降 雨 分 布 


Fig.3 Rainfall seasonal distribution in growing season in 2017 and 2018 


表 2 模型 参数 
Tab.2 The parameters of the threshold - delay model 


林 龄 /a k Ymax 7/d Rt max 
0~10 0. 384 0.312 7.442.8 0.2 0.107 
10 ~20 0. 836 0.403 8.3+3.1 0.2 0.236 
20 ~40 1.706 0.351 8.1+2.7 0.2 0.221 

>40 0. 892 0.349 8.2+3.2 0.2 0.117 


(图 4a) ,4 种 林 龄 梭 梭 NDVI 响应 速度 呈 递增 态势 ， 
其 中 林 龄 >40 a 的 NDVI 增长 率 变 幅 最 大 (10% ~ 
56% ) ,均值 为 30% ,反应 最 为 迅速 ;0 ~ 10 a 林 龄 增 
长 率 变 幅 最 小 (0.1% ~24% ) ,均值 为 12% ,反应 相 
对 迟缓 ; 当 降 雨量 为 2 ~5 mm 时 (图 4b),10~20 a 
的 变 幅 最 大 (11% ~ 83% ) ,均值 为 37% ,0 ~10 a 的 
变 幅 最 小 (12% ~24% ) ,均值 为 21% ,10 ~20 a 组 
NDVI 增 长 率 明显 快 于 0 ~10 a 组 ; 当 降 雨量 为 5 ~ 
10 mm 时 (图 4c) ,20 ~40 a 组 最 为 敏感 ,增长 幅度 
最 高 达到 170% ,均值 为 49% ,0 ~ 10 a 林 龄 变化 幅 
度 最 小 (2% ~ 37% ) ,均值 为 14% ; 当 降 雨量 大 于 
20 mm 时 (图 4d), >40 a 的 NDVI 变 幅 最 大 
(24% ~76% ) ,均值 为 57%,10 ~ 20 a 的 次 之 
(17% ~66% ) ,均值 为 36% ,20 ~40 a 的 增幅 最 小 
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(16% ~23% ) ,均值 为 18% 。 从 增长 率 变 幅 来 看 ， 
在 由 小 至 大 降雨 级 别 梯度 上 ,NDVI 最 大 变 幅 分 别 
为 20 ~40 a.10 ~20 a.20 ~40 a 和 >40 a 人 工 梭 梭 
林 ; 根 据 均值 则 分 别 以 >40 a、10 ~20 a.20~40 a Fi 
>40 a 为 最 大 。 这 表明 随 降雨 等 级 的 增加 ,NDVI 
增长 率 变化 最 大 的 林 分 不 尽 相 同 ,其 中 5 ~ 10 mm 
降雨 对 20 ~40 a 生 人 工 梭 梭 林 的 生长 发 育 影 响 最 为 
明显 ,0 ~2 mm 降雨 对 0~10 a 生 梭 梭 影响 最 弱 , 而 
林 龄 在 20 a 以 上 的 人 工 梭 梭 林 对 降雨 利用 能 力 更 
强 , 因 为 在 变 幅 和 均值 条 件 下 20 a 林 龄 以 上 级 别 出 
现 频率 达到 6 次 。 

从 不 同 林 龄 人 工 梭 梭 林 对 降雨 响应 时 间 来 看 
( 表 3、 图 5) ,降雨 发 生 后 植被 响应 不 是 立即 发 生 ， 
而 是 在 降雨 发 生 后 出 现 延 迟 ( 时 间 滞 后 )3 ~ 
15 d。 当 降雨 量 小 于 2 mm 时 ,大 于 40 a 的 林 分 响 
应 滞后 时 间 最 短 ,平均 为 7.2( +4.2)d;10 ~20 a 林 
龄 对 降雨 响应 时 间 最 长 ,响应 时 间 平 均 为 9. 3 
(+4.5) d。 当 降雨 量 介 于 2 ~5 mm 时 ,不 同 林 龄 
人 工 梭 梭 林 响 应 时 间 呈 递增 态势 ,0 ~ 10 a 组 响应 时 
间 平 均 为 6.3( +1.3)d,10 ~20 a 组 响应 时 间 平 均 
为 8.1( +2.1)d,20~40 a 组 响应 时 间 平 均 为 10.3 
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图 4。 不同 降 雨 等 级 下 NDVI 增长 率 随 林 龄 的 变化 特征 


Fig.4 Characteristics of NDVI growth rate with stand age under different rainfall level 


季 树 新 等 :不 同 林 龄 人 工 梭 检 林 对 降雨 的 滞后 响应 


表 3 降雨 事件 的 选择 及 相应 的 响应 时 间 
Tab.3 Selected rainfall events and the corresponding response duration 
响应 时 间 /d 
始 时 间 结束 时 间 降雨 量 /mm ul 
0~10a 10 ~20 a 20 ~40 a >40 a 
2017 -05 -13 2017 -05 -14 4.9 7 10 13 14 
2017 -06 -04 2017 -06 -06 27.1 8 5 7 6 
2017 -06 -20 2017 -06 -21 8.4 8 10 8 6 
2017 -07 -02 2017 -07 -03 9.2 12 11 8 11 
2017 -08 -01 2017 -08 -02 0.5 3 3 3 3 
2017 -08 -08 2017 -08 -09 3.8 5 6 7 7 
2017 -09 -12 2017 -09 -14 7.3 2 2 6 5 
2018 -05 -25 2018 -05 -26 1.2 15 16 15 13 
2018 -09 -04 2018 -09 -05 0.2 6 8 6 4 
2018 -09 -17 2018 -09 -18 6.1 5 8 4 
注 : 响 应 时 长 是 从 降雨 日 起 至 NDVI 达到 最 大 值 的 天 数 。 
18 p 00~10a x10~20a A20~40a o>40a 最 慢 ,响应 时 间 平 均 为 8.3( +2.8)d。 当 降雨 量 大 
15 ó o F 20 mm 时 ,10 ~ 20 a 组 响应 最 快 ,响应 时 间 平 均 
. pene 为 5.5 d;0 ~ 10 a 组 响应 最 慢 , 响 应 时 间 平 均 为 8 d， 
T ， x x 其 中 最 大 响应 时 长 为 16 d ,出 现在 2018 年 5 月 25 
£ i Q 、 ZA 土生 
Bi x p © g 日 降雨 之 后 。 比 较 4 种 林 龄 人 工 梭 梭 林 NDVI 对 降 
Eo. o 、 
ES 雨 事件 的 响应 过 程 可 以 看 出 (图 6) ,2017 年 6 月 20 
3} ð 日 发 生 8. 4 mm 降雨 ,不 同 林 龄 人 工 梭 梭 林 响 应 时 
0 i i i i i i 间 分 别 为 8 d(0~10 a) 10 d(10 ~20 a) .8 d(20 ~ 
-5 0 5 10 15 20 25 30 
降雨 量 /mm 40 a) 和 6 d( 大 于 40 a),0~10a 组 和 10~20a 组 
= i 为 逐 》 型 ,20 ~ 
图 5 降雨 量 与 响应 时 间 关 系 向 应 方式 为 逐渐 增加 型 ,20 ~ 40 a 组 和 大 于 40 a 组 


Fig.5 Relationship between RD of vegetation 


and rainfall events sizes 
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图 6 2017 年 6 月 20 日 降雨 后 不 同 林 龄 
RER NDVI 最 大 值 变化 


Fig.6 Maximum NDVI variation map of different ages’ 


Haloxylon ammodendron after the rainfall on June 20 ,2017 


(+3.3)d, >40 a 组 响应 时 间 平 均 为 10.9( +3.9) 
d。 当 降雨 量 介 于 5 ~ 10 mm 时 ,20 ~40 a 组 响应 最 
快 ,响应 时 间 平 均 为 6.8( +1.7)d;10 ~20 a 组 响应 


一 


向 应 方式 为 波动 增加 型 。 
3 讨论 


从 研究 时 段 代 表 性 来 看 ,研究 区 内 1981—2010 
年 多 年 平均 降雨 量 为 109.7 mm!" ,而 本 研究 选取 
时 段 是 2017 一 2018 年 ,年 平均 降雨 量 为 111. 95 
mm ,就 降雨 因素 而 言 ,本 研究 时 段 仍 处 在 平均 水 平 ， 
具有 代表 性 。 

从 降雨 量 在 主要 影响 梭 梭 NDVI 的 环境 因子 中 
的 地 位 来 看 ,不 同 林 龄 人 工 梭 梭 林 的 NDVI 增长 率 
与 降雨 量 、 气 温 ,干旱 间隔 期 (降雨 事件 前 无 雨天 
数 ) 的 Person 相关 性 分 析 表 明 ( 表 4) ,NDVI 与 降雨 
量 相 关系 数 分 别 为 0.546 .0.313 .0.345 和 0.338 ,说 
明和 干旱 区 人 工 权 权 林 的 NDVI 增 长 率 与 降雨 量 相关 
性 较 好 。 因 为 一 次 较 大 的 降雨 后 土壤 含水 量 增加 ， 
植物 光合 作用 增强 ,NDVI 增加 迅速 。 相 反 ,气温 与 
NDVI 增长 率 之 间 相关 性 差 , 其 诱发 原因 与 气温 升 
高 导致 植 被 叶片 气孔 关闭 、 光 合作 用 能 力 降 低 并 最 
AG Be NDVI FE BEA IE ,干旱 间隔 期 与 NDVI 增 


FOR 区 


10 ~20 a 


ie ee 
让 HD 十 
ae 


NDVI 增长 率 与 影响 因子 相关 关系 


Correlation between NDVI growth rate and its drivers at different habitats 


20 ~40 a >40 a 


表 4 
Tab.4 
0~10a 
降雨 量 0.546( P=0.102) 
干旱 间隔 0.347(P =0. 326) 
气温 —0.098(P=0.768) 


长 率 相关 性 较 差 ,虽然 干旱 区 土壤 有 一 定 的 “降雨 
记忆 ”, 但 生长 季 中 的 单 次 降雨 还 是 能 促使 植被 生 
产 力 在 较 短 时 间 内 快速 增加 ,其 中 0 ~ 10 a 组 NDVI 
增长 率 与 降雨 量 和 干旱 间隔 期 相关 性 较 强 ,这 是 因 
为 梭 梭 在 幼林 期 间 , 对 大 气 降雨 及 土壤 水 的 响应 更 
加 敏感 ,土壤 水 含量 越 高 ,生长 越 迅速 ("|。 

从 梭 梭 林 对 降雨 响应 期 来 看 ,不 同 林 龄 人 工 梭 
梭 林 对 降雨 的 响应 时 间 在 8 d 左右 (图 5) ,而 在 
2018 年 5 月 25 日 ,降雨 量 为 1.2 mm,4 组 梭 梭 的 响 
应 时 间 分 别 为 15 d,16 d、15 d、13 d, 说 明 在 生长 季 
早期 干旱 区 植物 生长 对 大 气 降雨 响应 并 不 敏感 ,这 
个 结果 与 戴 后 等 在 古 尔 班 通 古 特 沙漠 通过 主要 
建 群 种 梭 梭 生长 季 水 分 动态 研究 及 对 不 同时 期 相似 
量 级 (6.7 mm、7 mm 降雨 ) 降 雨 脉冲 的 响应 研究 观 
点 一 致 。 产 生 这 种 现象 的 原因 是 梭 梭 在 生长 季 初 期 
(6 月 之 前 ) ,水 分 需求 主要 依赖 于 土壤 中 冬季 的 解 
冻 水 ,此 时 由 于 气温 和 地 温 较 低 ,导致 梭 梭 在 生长 季 
前 期 生长 较 慢 。Li 等 "2 采用 时 间 分 辩 率 为 8 d 的 
MODIS 数据 ,探究 西北 和 干旱 区 沙 地 生境 和 荒漠 生境 
NDVI 对 单 次 降雨 的 响应 过 程 , 响 应 时 间 为 10 ~ 20 
d, 本 文采 用 时 间 分 别 率 为 1 d 的 MODIS 数据 , 较 之 
研究 结果 在 时 间 分 辨 率 上 有 明显 提高 。 降 十 发 生 时 
间 的 不 确定 性 ,使 得 梭 权 对 降雨 引起 的 土壤 含水 量 
变化 的 响应 也 具有 不 确定 性 ,进而 使 得 植被 NDVI 
的 响应 产生 不 同 的 变化 。 

从 不 同 林 龄 人 工 梭 梭 林 在 降雨 事件 后 的 响应 特 
点 来 看 ,在 经 历 2 ~5 mm 和 5 ~ 10 mm 等 级 的 降雨 
后 ,不 同 林 龄 人 工 梭 梭 林 NDVI 增加 的 波动 范围 明 
显 高 于 0~2 mm 和 大 于 20 mm 降雨 等 级 (图 4) ,由 
此 可 以 推断 2 ~5 mm 和 5 ~10 mm 降雨 对 人 工 权 梭 
林 的 生长 发 育 十 分 重要 。Sala 等 2 BE, PEAK 
地 区 <5 mm 降雨 对 植物 生存 和 生长 至 关 重 要 ,本 文 
认为 2 ~ 10 mm 降雨 能 促使 检 检 在 短期 内 迅速 萌 
发 ,促进 其 生长 和 NDVI 的 增加 。 不 同 林 龄 人 工 梭 
梭 随 林 龄 增加 高 度 、 冠 幅 都 呈 增 加 趋势 汪 ,虽然 梭 
梭 的 栽植 可 以 使 土壤 理化 性 质 发 生 一 定 程 度 的 改 


0. 313 (P =0. 378) 
0. 175 (P =0. 630) 
—0.087(P =0. 810) 


0. 345( P =0. 304) 
0. 023 ( P =0. 949) 
0. 334( P =0. 345 ) 


0. 338 ( P =0. 340) 
0. 187( P =0. 604) 
0. 020 (P =0. 956) 


变 , 但 随 林 龄 增加 土壤 有 机 质 含量 也 呈 下 降 趋 势 。 
本 文中 0 ~10 a 和 10 ~20 a 梭 梭 对 单 次 降雨 过 程 响 
应 时 间 相 较 于 其 他 2 组 更 短 , 可 能 与 土壤 有 机 质 速 
效 磷 等 含量 高 于 高 龄 权 梭 下 的 土壤 有 机 质 含量 有 
关 , 且 0~10 a Fil 10 ~ 20 a 梭 梭 根系 埋藏 深度 较 其 
他 2 组 更 浅 , 对 降雨 的 吸收 利用 更 早 。 

从 本 文 的 局 限 性 来 看 :@ 在 较 短 时 间 尺 度 内 
(2 a) 缺 乏 足够 的 有 效 降雨 ,导致 增幅 在 10% 以 下 
的 NDVI 动态 变化 过 程 往往 容易 被 忽略 , 且 研 究 区 4 
种 林 龄 人 工 梭 梭 林 之 间 在 植被 盖 度 、 密 度 及 存活 率 
等 方面 存在 差异 ,这 种 差异 在 现 有 的 研究 中 都 没有 
很 好 的 解决 。@ 本 文 研究 时 段 内 缺少 使 植被 NDVI 
增幅 明显 的 10 ~20 mm 及 大 于 30 mm 的 降雨 ,所 以 
对 这 两 种 级 别 降雨 响应 过 程 没有 体现 。@) 由 于 天 
气 原因 ,选取 独立 降雨 后 的 NDVI 连续 变化 过 程 中 
存在 有 云 和 云 影 的 情况 ( 出 现 无 效 数 据 ) ,导致 响应 
天 数 出 现 一 定 滞后 。 


4 结论 


在 不 同等 级 的 降水 条 件 下 ,不 同 林 龄 人 工 梭 梭 
林 表 现 出 不 同 的 响应 特征 。2 ~5 mm 和 5 ~ 10 mm 
降雨 是 影响 干旱 区 荒漠 - 绿洲 过 渡 带 不 同 林 龄 梭 梭 
生存 生长 的 两 个 重要 阔 值 范围 , NDVI 对 降雨 脉动 
的 响应 始 于 2 mm 左右 , 当 降 雨量 处 在 5 ~ 10 mm 
时 ,NDVI 迅速 增加 2% ~74% 。 降 雨 发生 后 ,不 同 
林 龄 人 工 梭 梭 林 的 植被 生产 力 响应 时 间 不 同 ,0 ~ 
10 a 组 为 7.4( +2.8)d,10 ~20 a 组 为 (8.3 +3.1) 
d,20 ~40 a 组 为 (8.1+2.7) d, 大 于 40 a 组 为 (8.2 
+3.2)d。 独 立 降 雨 是 引发 干旱 区 植被 生产 力 迅 速 
变化 的 重要 因素 。 
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Lag response of artificial Haloxylon ammodendron forests 


at different ages to rainfall 


JI Shu-xin, WANG Li-xiang, BAI Xue-lian, QIAO Rong-rong, CHANG Xue-li 


( Department of resources and environmental engineering ,Ludong University , Yantai 264025 , Shandong , China ) 


Abstract: Rainfall pulses are the most important constraints for the growth and succession of desert vegetation and 
directly affect the changes in vegetation composition, structure ,and function. As one of the dominant species of arti- 
ficial vegetation in the arid zones , Haloxylon ammodendron is a crucial factor in wind breaking, sand fixation , and 
water conservation. The objective of this study is to understand the response of the vegetation productivity of the arti- 
ficial Haloxylon ammodendron woodlands to rainfall events at different stand ages. We studied the Haloxylon ammo- 
dendron woodlands exhibiting different stand ages in a desert-oasis transition zone in the middle reaches of the 
Heihe River using the MODIS-NDVI (the approximate annual productivity) time series from May 1,2017 to Sep- 
tember 31,2017. Based on the corresponding daily rainfall data ,a threshold model and statistical methods were used 
to analyze the lag response process of the NDVI to different rainfall levels. The results showed that (1) When the 
rainfall was less than 2 mm,the NDVI growth rates in case of different stand ages were more than 5% and stands 
with ages of more than 40 year reacted the fastest; the growth rates fluctuated from 10% to 56% . When the rainfall 
was 2 -3 mm, the stands that exhibited ages of 10 -20 year showed the largest growth rate ,i. e. ,from approximate- 
ly 11% to 83%. When the rainfall was 5 - 10 mm, the most sensitive group had stand ages of 20 -40 year, with the 
highest growth rate being 170%. When the rainfall was more than 20 mm, the NDVI growth rate was greater for 
stand ages greater than 40 year than those for the remaining three stand ages, with an increase of up to 76%. (2) 
After the rainfall ,the NDVI response time was observed to be different for different stand ages: 7.4( +2.8) days 
for the 0 - 10 year group ,8.3( +3.1) days for the 10 -20 year group,8.1( +2.7) days for the 20 - 40 year 
group ,and 8.2( +3.2) days for the >40 year group. (3) Rainfall pulses contribute to rapid changes in vegetation 
productivity in the desert-oasis transition zone. Different stand age groups showed different response patterns to dif- 
ferent rainfall levels, and the NDVI increase for the artificial Haloxylon ammodendron woodland was the largest 
when the rainfall was 2 -10 mm. According to these results, when the desert-oasis transition zone is dominated by 
less than 10 mm of rainfall, the roles of 10 -20 year and 20 -40 year artificial Haloxylon ammodendron woodlands 
in the protection system are likely to be the most important factors. 

Key words: stand age; artificial Haloxylon ammodendron woodland; rainfall; NDVI; Threshold-delay model; 


lag response; Linze oasis 


